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I . 


Consider the ordinary differential equation 

X = f(t,x) (I) 

where x and f belong to R^, t is a real scalar and where f is assumed to 

be continuous on = {(t,x) :t^0, x <r},0<r^ + <». Let x(t,t^,x^) 

denote a solution of (I), such that x(t ,t ,x ) = x . In the following, if 

o o o o 

a, b both belong to and a^ b^ (a^ > b^) for all 1 = 1,2 nwe shall 

write a ^ b (a > b) . Further, a function f (x) where f and x belong to r’^, 
will be said to be of type K in a set S if f^(a) ^ f^(b) (1 = l,2,...,n) for 
any two points a, b in S with a^ = b^ and a^ ^ bj^ for all i = k. 

In addition to (I) we will also consider the "comparison" differential 
equation 

r = w(t,r) (II) 

where w and r belong to R™ and w is assumed continuous for all (t,r) such 
that t ^ 0 and r^ r”*. It will also be necessary to require w to be a 
function of type K, for each fixed t, in R™. 

The following lemma is given relative to (I) and (II) and is of funda- 
mental importance. It comprises the basic idea of the comparison principle. 

Lemma . Let v(t,x) belong to R™ and be continuous and locally Lipschitzian, 
with respect to x, on D^. Suppose that 

v*(t,x) = lim sup h ^ {v[t+h, x+hf(t,x)] - v(t,x)} ^w(t,v(t,x) 
h-K) + ~ 

for all (t,x) and also suppose v(0,x^) ^ r^. Let r^(t,0,r^) denote a 

maximal solution of (II) and let J denote its right maximal interval of 

existence. Then v(t,x(t,0,x )) < r (t,0,r ) for all t *J such that 

o — mo 

(t,x(t,0,x^))€ D^. 

The existence of the maximal solution r (t,0,r ), under the conditions 

m o 

placed on w, was established by Kamke. It can easily be shown that the local 



* 


2 > 


Lipschitz condition assures that v*(t,x) coincides with the upper right 
derivative of v(t,x(t,0,x^)) for every (t,x)6 D^. The lemma then follows 
immediately from another theorem of Kamke. With respect to the applications, 
it should be noted that 

V*(t.x) =1^ + 1^ - 

in the case v(t,x) has continuous partial derivatives. 

This lemma, or special cases of it, has been used by a number of indi- 
viduals to obtain sufficient conditions for the existence of solution pro- 
portion of classical stability theory. An example of its application can 
be found in the paper of Appendix I, "On a stability property of Krasovskil", 
which was published in the International Journal of Non-Linear Mechanics , 

Vol. 5, pages 507-512. In this paper, sufficient conditions for a type of 
practical stability, proposed by Krasovskil, were obtained. The results are 
of particular value to investigations involving the behavior of systems over 
a finite Interval of time. 

A different version of the comparison lemma was given in the paper of 
Appendix II, "A comparison lemma for higher order trajectory derivatives", 
published in the Proceedings of the American Mathematical Society , Vol. 27, 
Number 3, in March 1971. Here the basic lemma on comparison theory was used 
to obtain a comparison lemma which is useful when higher order trajectory 
derivatives of Liapunov functions are known. 

Appendix III, "A stability condition for linear comparison systems", 
published in the Quarterly of Applied Mathematics , July 1971, presents a 
result which is useful when applying the comparison principle following 
the basic approach given in a paper of Bailey. Using a result of this 
paper, it is possible to provide a very much simplified proof of an impor- 
tant section of Bailey's original paper. 



■ The paper of Appendix IV, ”A result on differential Inequalities and 
its application to higher order trajectory derivatives", was published in 
the SIAM Journal on Analysis , Vol. 4, Number 4, in November 1973. This 
paper was concerned with the important problem of obtaining a less restrictive 
condition for comparison systems than that imposed by the type K condition. 

In addition to obtaining a less restrictive condition, an example was given 
of its application to higher order trajectory derivatives. The example 
considerably improved that given in the paper of Appendix II. 

The result obtained in Appendix IV required investigation of the adjoint 
differential equation of the variational equation of the right side of the 
inequality. The condition itself could be considered as the failure to 
satisfy a certain two point boundary value problem. Consequently some effort 
was devoted to study of two point boundary value problems. Papers resulting 
from this investigation were those of Appendix V, "An existence theorem for 
boundary value problems", published in the Quarterly of Applied Mathematics , 
April 1973; Appendix VI, "Conditioning of linear boundary value problems", 
published in the Swedish Journal BIT , Vol. 12, 1972; Appendix VII, "Application 
of Liapunov theory to boundary value problems II", published in the Proceedings 
of the American Mathematical Society , Vol. 37, Number 1, January 1973; Appendix 
VIII, "Optimization of bounds for boundary sets", presented at the 1974 Inter- 
national Congress of Mathematicians , Vancouver, Canada. 

An application of the comparison principle to the analysis of approximate 
systems was given in the paper of Appendix IX, "Error estimates for approximate 
dynamic systems", to be published in a 1974 issue of the International Journal 
of Control . In this paper, two approaches were considered for obtaining 
estimates on the error between approximate and exact solutions of dynamic 




It sometimes occurs that the matrix of a linear comparison system is 
not of the required type K. In these cases it is sometimes possible to 
maintain the inequality while replacing the matrix by one whose off diagonal 
elements are the absolute values of the original matrix. The paper of 
Appendix X, "Stability of a class of matrices similar to nearly non-negative 
matrices", to be published in the SIAM Journal on Applied Mathematics in 
1975, is concerned with the question of when the stability of the former 
matrix implies the stability of the latter and whether a simple test for 
the stability of such matrices exists. 

Closely related to stability theory is the topic of bifurcation analysis. 
Roughly, a system will bifurcate when it loses stability. The papers of 
Appendices XI and XII are concerned with this problem. Appendix XI, "Sus- 
tained small oscillations in nonlinear control systems", submitted for publi- 
cation to the International Journal on Control concerns some results from the 
theory of bifurcating solutions which can be used to obtain conditions which 
allow sustained small oscillations in aircraft-spacecraft dynamics. Appendix 
XII, "Structurally stable bifurcating systems", submitted to the International 
Journal of Nonlinear Mechanics is concerned with showing that a system of 

m 

ordinary differential equations with a vertical bifurcation cannot be structur- 
ally stable. The purpose of these investigations was to determine the applica- 
tion of comparison theory to bifurcation analysis. So far, it is not clear. 

Finally, the paper of Appendix XIII, "Scalar n-th order comparison systems" 
submitted to the IEEE Transactions on Automatic Control discusses results from 
several papers and their application to comparison theorems involving higher 
order trajectory derivatives. 
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APPENDIX VII 


APPLICATION OF LIAPUNOV THEORY TO 
BOUITOARY VALUE PROBLEMS, II* 
by 

J. H. George** and R. J. York^ 

1. Introduction . In (1] , George and Sutton formulated a Liapunov 
theory for the existence and uniqueness of solutions to the following class 
of boundary value problems; 

(1.1a) x" = f(t,x,x*) 

(1.1b) x(a) “ A, x(b) * B 

where A,B are fixed n-vectors, ' *» d/dt, .x Is a n-vector with a 

continuous second derivative on [a,b], and the n-vector function f Is 

2n 

defined and continuous on D » [a,b] x R where (a,b] Is an Interval on 
the real line. The development In [1] replaced certain conditions of 
Hartman (2], by Liapunov- type results. Fundamental to this approach was 
the construction of a function g which Is equal to f on a bounded region. 

LEI'O^ 1.1. (Hartman [2, p. 432]). Supoose f Is defined and contin- 
uous on Djj ■ [a,b] x (x : | |x{ | <^ h} x r” and let N > 0 be given. Then 
for any E^x’^2^ there exists a continuous function g bounded on 

Dj^, such that f " g on •* (tj^,t2] (x ; | |xj | £ M) x (x* : j jx* | | £ N). 

By placing Liapunov conditions on 
(1.2) x'‘ - g(t,x,x’) 

It was shown In [1] that the solutions of (1.2) and (1.1b) are also 
solutions of (1.1). 

In this paper, a more restrictive class of Liapunov functions than 
those considered In [1] are Introduced. This allows the existence theorems 



la [1] to be proved with conditions In terms of f rather than g as In 

[1]. Thus, a principle difficulty In using the results of [1] are removed. 
Also, the new class of Liapunov functions Is still sufficiently large to 
Include Hartman's condition, [2]. 

2. Liapunov Conditions for existence of solutions . Let a Liapunov 
function V(t,x,x') be a continuously differentiable In (t,x,x*), real 
valued function.. Let 

3v .3v av. 3V ^ V 

3x " ^3x ’•••*3x^‘ 3x* “ ^3x;’‘*‘»3x’' 

In In 


Let <• ,«> denote the dot product In n-dlmenslonal Euclidean space. Then 
define the derivative of V with respect to t along a solution of (1.1a) 
to be. 


2^!!. 4. v’> + < 

at ^ 3x»* ^ 
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Also let x(t) belong to the class C (tj^,t 2 l If x’ (t) is continuous 
on [t^,t 2 ]. We will nov; restate a Lemma of George and Sutton [1] In terms 
of this Liapunov function. 

LEMMA 2.1. (George and Sutton [!)). Suppose there exists a Liapunov 
function U(t,x,x') defined on D 2 = [tj^,t 2 J ^ (x: 1 1 x| | £ m} x r”. Let 
X(t) be the set of functions x(t) e C [tj^,t 2 l satisfying ||x(t)|l 
on [tj^,t 2 l, and also the following properties: 

(1) U(tj^,x(tj^) ,x’(tj^)) -0 

(11) U(t,x,x') ^ (t - tj^) (||x'H), where 4> Is a positive continuous 
function defined on [0,®») such that <|>(r) ■► “» as r 

(Ill) there exists a positive constant L such that 

U' = ^ interior of D 2 . 



? 


Then there exists a positive constant N such that ||x'(t)|| ^ N on 
Iti.t^l, for every x(t) belonging to X(t). 

EXAMPLE. Let 

U - (t - tj^)j|x’||^. 

Then conditions (1) and (11) are satisfied. Condition (111) becomes 


(2.1) <x'(t),x"(t)> ~2(i''- "t^) ^ 


on 


(t » ^2 1 


From the lemma, there exists an N such that If (2.1) Is satisfied and 
|}x(t}|{ then ||x'(t)(| jc N on [t^^.t^J. Many lemma’s similar to 

Lemma 2.1 could be developed. This was simply an attempt In [1] to show 
that alternatives existed to the standard Nagumo type theorems [2, p. 428-9]. 

THEOREM 2.1. Assume there exists a Liapunov function satisfying the 
hypothesis of Lemma 2.1. Let V(t,x,x’) be another Liapunov function 
defined on " (tj^,t 2 ) (x: |{x|| ^ M) x r” and satisfying: 

(1) <1^, ,f><0 for ||x||>M, j|x’||^N+l 

and 


Then vl > 0 In the Interior of D, Implies v' > 0 In the Interior of 
f — 3 g — 

Proof. Let N be the constant provided by Lemma 2.1. As In Hartman 

I 1 - 6(jjx’j| - U)f(t,x,x’) on 


(2], let g(t,x,x') 



Here, 


6(s) 


i“ 


If 

If 

if 


x < 0 

0 _< 8 j< 1 

s > 1 


The proof now consists of establishing that 


V > 
8 - 


0 on D with 


m 


= 11 *= M. 



Then, the argument Is extended to ||x|| > M. The proof is motivated by a 
result of Knoblock [3]. 


On Djj with I |x{ I - M. 


( 2 . 1 ) 


3V , ^3V . . ^ 


where t ■ 1 - £(||x'|| - M) e (0,1J. Since (2.1) is a linear function of 

T £ [0,1], if V >0 for T = 0 and T ® 1 then V > 0 for all 
8 — 8 — 

T £ [0,1] . At T " 1, Vg « Vj ^ 0. At T ■ 0, | |x' 1 1 > N + 1, we have 

from condition (i) 




If I|x|| > M, 


„• _ 3V . .3V . ,av M Ifac 

\ ^ ^'*TRr 


9V . ,3V . !IT ,3V Hx . „ 

3t * ” * rrar 'k- •‘'‘’Trar-’' ^ - ® 


by repeating the argument above for T «» 0 and T = 1, using condition (11) 
when T « 1 and condition (1) when T ■ 0. 

Using this result, THEOREtl 4.2 in [1] can now be stated as follows: 
THEOREM 2.2. Let x(t) be a solution of (1.2) defined on [t^^,t 2 ] 
[a,b] and suppose that ||x(t^)|| M, ||x(t 2 )I| Let there exist a 

Liapunov function V(t,x,x') defined on such that 

(i) V(t,x,x’) ® 0 whenever j|xjj » M 
(ii) V(t,x,x*) > 0 whenever Ijxj) > H 

(ill) < 1 ^ 1 . f> 1 0 foe lUn > IU llx’JI i K + 1 

<*''> for (|x||>M, ||x'|li» 

(v) > 0 in the interior of 



Then ||x(t)i| £M for all t e 

t 

Proof. From Theorem 2.1, condition (v) can be replaced by V ^ 0 In 

8 

the Interior of D^. The proof continues exactly the same as the proof of 
Theorem 4.2 in [ 1 ] . 

2 2 2 

EXAIlPLE. V = Ijx'll (||x|| - ) shows that need not be zero. 

In this case, (ill), (Iv) and (v) will impose conditions on f. 

THEOREII 2.3. Suppose the Liapunov function in Lemma 2.1 has the 
additional properties ^ L In the Interior of D 2 and 

(1) <|^,,f>>0 for |jx|(<Il. l|x’|llN + l. 

t 

Then U < L In the Interior of D-. 
g - 2 

Proof. In a similar fashion to Theorem 2.1, 


Ul - <||,x»> + <|^..f> T 


g 3t 

fl 


~ 9t 


‘3x> 


At T - 1, - u’ < 0. When T = 0, Up = < |r + + 

dll t 

<^, ,f> " Uj L, by condition (1). 


g 3t 


The results of [1] can now be stated entirely in terms of f. 

THEOREM 2.4. Suppose f(t,x,x') Is defined and continuous on 
2x1 

[a,bl X R . Suppose there exists two Liapunov functions V(t,x,x') and 
U(t,x,x') as described In Theorems 2.2 and 2.3. Then the boundary value 
problem (1.1) and x(tj^) “ x^^, x(t 2 > * * 2 , ||xj^||, IIX 2 II has at least 

one solution. 

Proof. From Theorem 2.2 and 2.3, we have v' > 0 in the Interior of 

8 — 
t 

Dj and U^ jc L in the Interior of D 2 . The proof is now the same as 
Theorem 4.4 In [Ij. 

EXM'IPLE. Hartman [2, p. 433] gives the following condition to Insure 
||x(t)|(_<ll on [tj^,t 2 ]: 



<x.f> + I |x' 1 1^ > 0 If <x,x’> “ 0 and | |x| | ^ M 

By choosing V(t,x,x') « <x,x> - * 2<x,x*>, V^' = 2[<x,f> + 

2 

jjx'J) ]. Hartman's condition Implies V evaluated along a solution x(t) 

of (1.1a) does not have a maximum at any point t c [tj^,t2l when 

||x(t)|| ^ It* This follows since If “0, > 0* I^ow V = 0 when 

I |x| I » M and V > 0 If | |x| { > M If there existed a solution x(t) 

such that I |x(t) 1 1 > M, then ^ 0 and all conditions of the theorem 

3V 

are satisfied with the observation that ° 0 
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APPENDIX VIII 


VALUE PROBLEMS* 

1. Consider the boundary value problem given by 

(1) y" = f(t.y.y’) 

(2) y(0) = A y(l) = B 

2 

.where f(t,y,y' ) is continuous on [0,1] X R and A and B are real. 

A recently developed technique for solving boundary value problems 

with very general endpoint constraints developed by Bebernes and Wilhelmsen 

[IJ* [2], [3], did not lend itself directly to boundary value problems 

(BVP) (1) “ (2) . Beyond the assumption of lower and upper solutions of 

equation (1) their principal restriction was a differential boundedness 

condition imposed on the solutions of equation (1) . 

By strengthening the differential boundedness to a Nagumo condition 

Bebernes and Fraker [4] demonstrated a modification of the techniques 

suitable for solving BVP (1) - (2) . The results possible because of the 

modification have been discussed in different settings by Jackson and 

Schrader [5] and Jackson and Klaasen [6]. 

Examples are given showing that the differential boundedness condition 

is neither sufficient nor necessary for demonstrating the existence of 

solutions to BVP (1) - (2) using the techniques in [1], [2], or [3]. 

2 

2. Let a(t) and 3(t) e C [0,1] be lower and upper solutions, respectively, 
of equation (1) on [0,1] with a(t) < 3(t) for ail t e [0,1], A solution, 
y(t), of equation (1) is said to be initially differentially bounded 


EXAMPLES CONCERNING SOLUTIONS TO BOUNDARY 

Ross Fraker 



(to the right) with respect to a lower solution - upper solution pair 

(a < 0) if there exists a positive N(n,tQ) such that. |y'(t)( < N(n,tQ) 

for t e [0,tQ] whenever |y*(0)( < n and a(t) < y(t) < g(t) for all t e [O.t^]. 

Let be a compact connected set in the (0,y,y') - plane such that 
there exist points (O.y^^.yj^') and (0,y2»y2') with y^^ = 0(0), y^' > 0*(O), 
y 2 = o(0) , and y 2 * < a* (0) , and let S 2 be a closed connected set in the 
(tQ»y»y') “ plane such that the set 

S 2 n {(tQ,y,y') : “(t^) < y < 3(tQ>, y' = k (arbitrary)} 

is non-empty for every k. If the solutions of equation (1) are initially 
differentially bounded with respect to some lower solution - upper solution 
pair (a < 0) then there exists at least one solution to equation (1) with 
(0,y(0),y' (0)) e (t^.y (t^) ,y * (t^)) e S 2 , and o(t) < y(t) < 0(t) for 

all t e [Oyt^]; (see the theorem of [3]), It is clear that this result 
can not in general be applied to BVP (1) - (2) . 

3. It is noted in [2] that the properties of the set are not necessary 

for the existence of a solution originating between o(0) and 0(0). It 
is sufficient only that the solution funnel of a compact connected component 
of the initial set Intersect both the upper solution and the lower solution 
at least by time t^. A Nagumo condition can be used to eliminate the neces- 
sity of checking this property. Under both conditions BVP (1) - (2) can be 
solved. The following example shows that even the weaker differentially 
boundedness condition is not necessary for the solvability of BVP (1) - (2). 

Consider y" = (y’)^ on [0,1]. Let a(t) = -1 and 0(t) = 1. In order 

to solve BVP (1) - (2) for a(0) < A < 0(0) and a(l) ^ B < 0(1) it is suffi- 
cient to take = {(0,A,y’) ; (y'( < 2} since by the mean value theorem 



/Q> 


y'(t) > y'(0) + t-min f(t,y,y') > 2 
0<t<l 
-l<y<l 
y*>2 

when y' (0) > 2 for example. The connectedness of the funnel cross section 
at t B 1 yields a solution If 

Sj = {(l.B.y*) : |y'| < «»1 . 

The example we give now shows that differential boundedness is not 

sufficient to solve BVP (1) - (2). 

3 

Consider y" = -(y*) on [0,1]. Let o(t) = -1 and 6(t) = 1. Clearly 
the solutions, when they exists are differentially bounded. Let A = -1 
and B = 1. To solve the differential equation we consider the auxiliary 

3 

equation z' = -z to find that 

y’(t) = (2t + (y'(0)) , y'(0) > 0, 

• 1/2 

y^t) a -(2t + (y'(0))”^) , y'(0) < 0, and 

y'(t) =0 , y'(0) a 0. 

Since we want y(0) = -1, y(l) = 1, and -1 < y(t) < 1 for all t e [0,1] we 
must take y'(0) > 0 and because of the form of the equation we must require 
y' (0) > 2. This leads to the solution 

-2 

y(t) = (2t + (y'(0)). ) +c 
for appropriate constants y* (0) and c. 

The boundary values lead to the equations 

-1 a |y'(0)| ^ + c and 
_9 1/2 

1 a (2 + (y' (0)) + c. 



Eliminating c we find that 

1 9 1/2 

2 + |y’(0)r-^ = (2 + (y'(O))"'^) 

or 

4 + (y'(0))"^ + 4|y'(0)l = 2 + (y'(0))"^. 

Hence we must have 

|y'(o)|=-| 

which is clearly impossible. Thus no solution to this boundary value 
problem is possible. 



APPENDIX I 

STABILITY OF A CLASS OF MATRICES 
SIMILAR TO NEARLY NON-NEGATIVE MATRICES 

R. W. Gunderson and E. E. Underwood* 

Abstract . A simple stability condition is obtained for a class 
of matrices which occur in the application of linear differential inequalities 

n ^ 

1. Introduction . Let a,b belong to R and write a b if 
and only if a^ ^ b^ for each i = l,2,....,n . Let A be an n x n 
real matrix and z(t) a vector function satisfying 

(1) z' (t) Az(t) 

for te[t^,T]. Suppose further that 

z(t ) < x(t ) 
o — o 

where x(t) is a solution of. 

x' = Ax 

over the same interval [t ,T]. Then it is known ([l],p.l37) that 

o 

z(t) ^ x(t) 

for all tE[t^,T] if and only if the matrix A is nearly non-negative; 

that is, a^j ^ 0 for all i j and i,j = 1,2 n. Further, it is 

known that the nearly non-negative matrix A is stable if and only if its 
k^^ successive principal minor has sign (-1)^ ([2],p.-74). 

The above results have a number of interesting applications, ranging 
from estimating the errors of approximate solutions of differential equations 
[3] to the stability analysis of composite, or inter-connected, 

. *Dept, of Mathematics, Utah State University. The research of the 
first author was supported by NASA Grant No. NGR-002-016. 
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systems [A]-[6]. It often occurs, however, that the matrix A of (1) 
does not turn out to be of the required type, i.e. nearly non-negative. 

In these cases it is sometimes possible to maintain the inequality (1) 

and at the same time replace A by |aJ , where |a1 is obtained from A 
by replacing the off-diagonal elements by their absolute value. Thus, 

(a] is nearly non-negative and the preceding results apply. 

It is of Interest in the applications to know when the spectrum of 
A is retained by 1 a|. In particular, it is of interest to know when the 
stability of A implies the stability of 1 a 1 and whether a simple test 
for the stability of such matrices exists, such as it does for |a| . 

2. A Family of Matrices Similar to Nearly Non-Negative Matrices . 

For a given n x n matrix A, draw a line through the first row in which 

a negative non-diagonal element appears. Denote the number of the row by 
ij^ and draw a line through the column i^ . Now draw a line through the 
next row which has a negative non-diagonal element, not in colimin ij^. 

Denote this row by I 2 and draw a line through column i^, and so forth. 

In this manner the matrix is covered with lines over every negative element. 
Now suppose that every non-diagonal element covered by a line is non- 
positive, with the exception of those covered at intersections of lines, 
which are non-negative. Then it is clear that there exists a permutation 
matrix P such that 

(3) /A B 

P AP = { 

\®2 ^2 



:re the square matrices A, and A» are nearly non-negative and 



are non-positive. If is of dimension (n-m) x (n-m) and A 2 is 

m X m, then A will be similar to a nearly non-negative matrix K, where 



Remark 1 . In the applications of comparison theory, (1) occurs as 
a system of differential inequalities satisfied by the trajectory 
derivatives of vector Liapunov functions [3]-[6]. Since the component 
Liapunov functions are almost always postive definite, it is worthwhile 
to consider (1) under the assumption z ^ 0 and with A of (1) assumed 
similar to a nearly non-negative matrix as described above. In such cases 
define the vector y by 



Since y is only a rearranged version of z , it follows from (1) and 
the definition of P that 


< p'^APy 



and, since y ^ 0» finally that 


l.e. 


y' 1 Ky 
z' < IaIz 


where 



Remark 2 . In equilibrium theory of economics, a matrix A satisfying 
(4) is called a Moroshima matrix if, in addition, it is irreducible and 
has all negative diagonal elements [7]. The requirement of irreduclbility 
allows a generalized Perron-Frobenius theorem to hold, from which some 
of the results of this paper can be deduced. The approach taken here does 
not require the Perron-Frobenius theorem and, hence eliminates the restiction 
to irreducible matrices. 

3 . Stability of Matrices Similar to Nearly Non-Hegative Matrices . 

Following Marcus and Mine [8], let ^ denote the set of all 

strictly increasing sequences of length r, formed from the first n: 
integers. Any submatrix of the n x n matrix A can then be represented 
by A[a|8] where a and 3 both belong to Q for 1 ^ n and 

where a determines the sequence of rows and 3 the sequence of columns. 

Using this notation, an n x n matrix H is said to be Hicksian [7] if 

det A[a|a] = (-1)*’ jdet A[a|a] | 

for each aeQ and all 1 ^ ^ ^ ® > that is, if every r^^ order principal 
r , n 

IT 

minor of the matrix A has sign (-1) . 

Definition . A n x n matrix P will be called a signed permutation 
matrix if the matrix obtained by replacing the elements of P by their 
absolute values is a permutation matrix. 

Lemma . Let P be an n ^ n signed permutation matrix and 
consider the family of all r x r minors 

det P[a)3] 

a,3eQ . For a given a there corresponds one sequence 3* 

I* «n 


where 


such that 


det P[a|3*] = ± 1 


and 

det P[o|b] *= 0 


otherwise. 

Proof . First observe that 

PP*^[ala] = 1 

and that 

det P*'^[w|a] = det P[a|w] 


for any “,weQ . By the Cauchy-Benet theorem 


r.n 


det PP^[a|a] = ^ det P[o|w] det P^[w|a] 

weQ 


'r.n 


I 

WeQ 

^r.n 


(det P[a|w]) = 1 


It follows that det P[a|w] “ +1 for one and only one w and is zero 
otherwise. 

Theorem 1 . Let A be Hicksian and P a signed permutation matrix 
both of dimension n x n. Then, the matrix 

H = PAP'*^ 


is also Hicksian. 



Proof . Again apply the Cauchy-Benet theorem to obtain 
det PAP^[a]a] =» det PA[a|w] det P^[w]a] 

weQ„ „ 

r.n 


^ det p"^[w|a]^, det P[a|y] det A[y|w] 


weQ YeQ 

r>n ^ ^r,n 


weQ yeQ 
r,n ' r,n 


det P[a|w] det P[a|Y] det A[y|w] 


By the lemma. 


det P[alw] det P[o|y] “ 0 


unless w = Y = And It follows that 

det PAP^[a|al = det A[8*(3*]. 

T 

Since A is Hicksian, so is PA? . 

Theorem 2 . Suppose there exists a signed permutation matrix P 
such that 

PAP^ = K 


where K is nearly non-negative. Then |a| is stable if and only if A 
is Hicksian. 

Proof . From the similarity established in the preceding section, it 

follows that if |a| is stable then K is also stable. Since K is stable 

k 

its successive principal diagonal minors have sign (-1) , as mentioned 
in the introduction. By a result of Kotelyanskil ([2],p.71) it follows that 
K is Hicksian. Hence by the lemma, A is Hicksian. If A is Hicksian, 
then K is Hichsiah and thus stable, since its successive principal 
diagonal minors then have sign (-1) . But, if K is stable, |A| is stable. 
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1. Introduction 


Consider. the system of first order equations 

X* = f(x,u) (1) 

where y is a parameter, f is assumed analytic in x and y at the origin 
of (x,y) space and f(0,y) = 0. Suppose there exists a family of 
solutions 

x(t,y) (2) 

such that to each neighborhood of the origin (x,y) there corresponds 
at least one value, y, such that (x(t,y), y) is contained in that 
neighborhood. In this case, the origin is called a bifurcation point, 
the family of solutions x(t,y) is called a bifurcating branch and the 
solutions corresponding to fixed values of the parameter y are called 
bifurcating solutions (Krasnosel'skii 1964). 

Example. Consider the scalar equation 

x' = x(x - y) 

It can easily be seen that 

x(t,y) = y 

is a bifurcating branch of solutions corresponding to a bifurcation 
branch at -the origin of the two-dimensional (x,y) space and that the 
bifurcating solutions are constant, or steady state solutions. 

In the following we shall utilize some results from bifurcation 
theory to investigate the existence of small amplitude periodic behavior 
in launch vehicle dynamics. It will be assumed that the nonlinearity 
exists as a cubic term in the rudder response. 
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Starting with Poincare, there have been , a number of important 
contributors to the theory. Among the early contributors to the 
theory of periodic bifurcations were Hopf (1942) and Friedrichs (1965) . 

We shall follow quite closely the approach given in Sattinger (1973) . 

2. Bifurcations in System Theory 

In addition to the existence of bifurcating solutions, either steady 
state or periodic, it is usually necessary in practice to determine their 
stability properties. The definition of asymptotic stability of steady 
state solutions is well-known and need not be presented here. However, 
the definition of orbital stability of periodic solutions is perhaps 
less well-known and given here for convenience. 

Definition . Let Y denote the closed path x = p(t) in x-space. 

The periodic solution p(t) is said to be orbitally stable if for each 
e > 0 there exists a corresponding 6 > 0 such that every solution 
x(t) of (33) whose distance from y is less than 6 for t = 0 is 
defined and remains at a distance less than E from y for sll 
t ^ 0. It is said to be orbitally asymptotically stable (and y is 
said to be a limit cycle) if in addition the distance of x(t) from y 
tends to zero as t 

Consider now the special case of equation (1) given by 

x’ = (A + viB)x + N(x) N(0) = 0 (3) 

where A,B are real n x n matrices and N(x) is an analytic vector 
function beginning x^ith non-zero k-th order terms, k 2. In addition, 
assume that the pair (A,B) has one of the following properties: 



Property I . The matrix A has a simple eigenvalue at the crigin and 
T 

if A(J) = A ip =0, then <B4> > r 0, where <*,*> is the dot product 

o o o o 

of r’^. Further, the matrix A + yB has all roots with negative real 

part when y < 0 and has exactly one unstable root for y > 0. 

Property II . The matrix A has a pair of simple, pure imaginary, roots 

T 

at ±i(jd and if A? = iw ,A ri = iuri then <B£ ,ri > 0. Further 

o o o o o o o 

assume A + yB has all its eigenvalues with negative real parts when 
y < 0 and has exactly one unstable complex pair if y > 0. 

The following theorem can then be given relative to the above 
properties: 

Theorem . (Hopf (1942), Sattinger (1973), Friedrichs (1965)). 

a) Assume that (3) satisfies Property I and, in addition, assume that 

<N(<{) ),ij; > ^ 0. Then the origin is a bifurcation point and the 
o o 

bifurcating branch consists of asymptotically stable steady state solutions 

b) Assume that (3) satisfies Property II and, in addition, that 

<N(C ),n > # 0. Then the origin is a bifurcation point and the 
o o 

bifurcation branch consists of orbitally asymptotically stable periodic ' 
solutions. 

Remark. The theorem is proven by showing the existence of an analytic 

one-parameter family of solutions (x(e),y(e)) satisfying x(0) = 0, 

y(0) = 0. The conditions <N(4> ),<p > 0 and <N(C ),n > 0 then 

o o o o 

allow the use of the implicit function theorem to solve for e as a 
function of y, providing the bifurcating branch. 

Re~:ark . Properties I and II both require that the simple eigenvalues, 
either real or pure imaginary, pass from the left to the right side of 



the complex plane when y = 0. This, of course will, not occur in general, 
but a simple transformation will place the system in proper form. 

Suppose, for example, the crossing over occurs for y = y^ in equation • 
(3) . Then set 

X = y - y 

o 

At = A + y B 
1 o 



and the new system 

x' = (A^ + XB^)x + N(x) 
will satisfy the condition. 

3 . Application to Flexible Body Dynamics 

Consider the following system of equations, obtained from the analysis 
of flexible body dynamics together with rigid body and control system 
dynamics . 

^r + ^ 

•• *22 ' 

<i), + 2k.?W9. + 0) (|), = to k 3 

bob 

3 

6 = f(a) = a - a^cr 
cf = + <|)^) + 

where 

(|)^ - attitude angle of the rigid body 
3 - rudder (control force) deflection 

C 2 “ control effectiveness coefficient 
a^, ~ control system gains 

- attitude angle due to bending dynamics (fixed mode) 
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5,oj - damping and natural frequency of bending mode 
k - normalization constant 
a - control command 
f(0) - nonlinear control command 
Letting 



equations (4) can be put into the form 

x' = (A + yB)x + ba^ 

where 




3 / 


and 

y = +k 

It can easily be seen that there exists a pyre imaginary pair of complex 

eigenvalues at rjw , providing *^2^1 ^ remains to enforce the 

conditions of Property II. To find ^ and r) it is necessary to 

o o 

solve 

“ ±o)E 
o o 

A% = iwri 
o o 

Carrying out the calculation, it turns out that 


A 




/“\ 


’^o “ 


0 

1 

ioj 


I 


where b^, b^ are complex constants depending on a^, w and C 2 » 
Consequently, the dot product conditions become 


and 


if 


<B^ ,n > = 2 Cw 0 

o o 


<N(C ),n > / 0 

o o 


B2 ¥ 0 

3a?w^ (a b- + a + a^b^)^ 
1 ' o 1 0 12 

a b^ + a + a.b„ ^ 0 
o 1 o 12 


4. Conclusion 


In the preceding, we have assumed a cubic nonlinearity in the rudder 
dynamics and have determined conditions under which a bifurcating branch 



of orbitally stable periodic solutions will exist. In the case 
considered, it was possible to determine rather easily conditions under 
which the system matrix had a pair of simple, pure imaginary, eigenvalues 
In more complicated cases this can still be accomplished by utilizing 
various linear stability techniques. The D-decomposition method of 
determining stability regions (Siljak 1969) ought to prove especially 
useful in this application. 
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1. Introduction . In the applications it is usaully required that a 
system possesses the property of structural stability. Briefly, this 
means that small changes in the differential equation result in 
topologically equivalent solution behavior. In this paper we show that 
a system which has a vertical bifurcation cannot be structurally stable 
in a neighborhood of the origin and develop a procedure to determine 
the existence of such bifurcations. 

2. Structural Stability and Vertical Bifurcations . 

Consider the system 

X* = <A+Xb)x + N(x) (1) 

where A,B are n x n matrices and N(x) is an analytic vector 
function defined or. and begins with non-zero k-th order terms, 

k >. 2. We shall . iume, for simplicity, that N is a bilinear form. 

Suppose there exists a family of .soluti >ns 

x(t,X) I 0 

such that to each neiohLcrhooc of the origin of the n + 1 space 
(x,X) there correspouac: at least one value, X, such that ^(t,X),xJ is 
contaiuec .n that .reigh jorhood. Then the origin is called a bifurcation 
point; the family of solutions x(t,X) is called a bifurcating branch 
and the solutions corresponding to fixed values of the parameter are 
called .if urcating solutions (Krasnosel'skii [1]). 

A system will be called structurally stable in a neighborhood 
of the origin if, for any sufficiently small perturbations, there exists 
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a hotneoivorpKism carries each trajectory of the original system to 

a trajectory of the perturbed system (Hirsch and Smale [2j). A 
bifurcating branch (x(t,A),X) is structurally stable if for all values 
of X 0 , sufficiently small, the system obtained by substituting 

X = y - x(t,X) into (1) is structurally stable in a neighborhood of 
the ^ L _^gr A4 • 

3. gxletence o-f Vertical Bifurcations . Folloxving Sattinger [3,4], 

assume (i) that A has a one-dimensional null space, Au^ = 0. Since A 

is n >r n matrix it satisfies a Fredholm alternative. That is, the 

* 

system Au = f has a solution if and only if < > = 0, where 

A u =0 ani < u _,u > = 1. Further < Av,u > = 0 for all v in 

" o 0 0 o 

r I n * I 

Ie_ M = tu u G R , < u,u > = 0/. 3y the above A J.S invertible 
o o 

on M , Denote its generalized inverse by K an<i the projection onto 
o o 

by Q. 

Finally, assume (ii) A + XB has its spectrum in the left half 
plane for X <0, h . a one dimensional null space when X = 0 and, 
when X < 0, this tigenvalu.e ‘crosses into the right half plane, while 
all others remain in the left half plane. 

o 

The steady state bifurcations of (1) can now be obtained by 
setting 

(A + XB)u + N(u,u) = 0 (2) 

c * 

* 

Let CyJ denote the inner product < ^ and introduce a parameter 

e into (2) by setting 



where 


C = Q(u^+0. [?] = 0. 

Equation (2) then takes the form 

6 

AC + eaB(u^+C) + eN(u^+C,u^+C) =0. (3) 

Operating on (3) first by KQ, then taking the [•] on (3) leads to 

C H- eoKQB Cu^+C) + eKQN(u^+C,u^+C) =■ 0 (4) 

and 

a[B <J)^] + cf[B Cl + [N(u^+C,u^+0] = 0 (5) 

Equations (4) and (5) are called the Liapanov-Schmidt equations [3,4]. 

The implicit function theorem can now be used to prove the following: 
Theorem 1 . (Sattinger [3,4]). Under the assumptions (i) and (ii) , and 
if [B u^] 0 there exists an analytic one-parameter family u « u(e), 

X = X(e) of solutions of (2) such that [u(e)J = e and X(0) = 0, 
u(0) =0. 

Remark 1 . The statement of the theorem and the proof do not rule out 
the possibility that u(eX f 0, X(e) = 0. If this is the case, the 
corresponding branch is cttller<l a vertical bifurcation. 

Exam ple Ctmsicter the scalar linear equation 

X* ■ Xx 

With each X = 0, the equation for steady state solutions is satisfied 
for any value of x. That is, a vertical bifurcation occurs at the 
origin, which is a bifurcation point. 
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which also satisfies the conditions of Theorem 1. The steady state 

2 

solutions are u = 0 when X 0 and U 2 = when 

X =0 and / 0 . 

•We next show that a vertical bifurcation cannot be structurally 
stable. 

Theorem 2 . Suppose the system 

X* = F(x) (7) 

has the property that 

(a) there exists G^(x^) such that the solution x^ = 0 of 

= F(Xj^) + G^(x^) is asymptotically stable 

(b) there exists G 2 (x 2 ) such that the solution X 2 = 0 of 

X 2 = I^(x 2 ) + G(x 2 ) is unstable. 

Then the system (7) is not structurally stable. Here F, G^, G 2 are 
sufficiently smooth, 1 g^(x) 1 < a |x|, and a sufficiently small. 

Proof : Suppose the system is structurally stable. ■ It follows that their 

exists a homeomorphism X 2 = Tx^^. Since the system of part (a) is 
asymptotically stable, its solutions form a compact set so that the solutions 
of part (b) must form a bounded set. However, the solutions of (b) 
are unstable and hence not bounded. 

Corollary . If the system (2) has a vertical bifurcation then it is 
not structurally stable. 



Exigtonce of Vertical BifurcatJ-ons . Assume that 

2 

C = * 

a = Oq + + ••• 


and substitute into equation (4) to obtain 
= kqn((!)q.4)q) + OqKQB (j>Q 

-^2 = 2KQN(4)q,?^) + a^KQB (()q + OqKQB 




c. 


-^n “ 2 KQN(<?q,C^_P + ••• + 


/2KQN(C^^ , , n 

2"^ 1 

KQN(C ,, ,), n 

. n-i r.-^x 


even 


odd 


(S) 


"n-lP ••• 


anc .. o (5) to obtain.the solvability conditions: 

(Pq] + [N(4>q,4)q)] = 0 

0^[B <J»q] + 0q[B Cj^] + 2[N((j)Q,q)] = 0 

: ( 9 ) 


Vl'“ *oI 






,^2[N(^^ , C^)] , n evei?j 

i ' ,!■ 

I ?^_l)], n odd j 

2 ~2~ 


= 0 



"Page missing from avaiiabie version" 



Theorea 2 . Suppose [N(c,,(J;„)] = 0 and 

■ <J U 

and that there exists a basis {x^.x^, .. 

space spanned by ••• > ...) 

1 ^ n 

[N((})q,x^)] =0 i = 1,2, . 

[N(x^,x J ] = 0 i,j = 1,2, 


is not 

u u 

• , x^ i, ic^n*“l, 

ic — 

which satisfies 


. . ,k 


i < k 


identically zero 
for the 


( 10 ) 


Then equation (8) and (9) have the solution 5(t;) y 0, 0(e) = 0. Tliat 
IS, (2) has a vertical bifurcation. 

Proof : Since {x^, ... is a basis and since N(x) is a bilinear 

form it follows that condition (10) implies for i > j, d,j = 1,2, ••• 


■N((^Q,q)] = 0 

CN(Ci,C.)l = 0 


(li) 


From (9), o(e) = 0. Solving (8) for ’ yields -the desired 
nonzero solution. 

Theorem 3 . There always exists a family of bilinear forms for which the 
hypothesis of Theorem 2 is satisfied, providing [N((f>Q,cj)Q)] = 0 and 
N(4)Q,yQ) is not entically zero, when n ^ 2. 

Proof; From (9) = 0 and from (8) 

V 

= -KQN((i)Q,(j)Q) . 


Let 


and assume 



Then, from ( 8 ) 


= -2KQN(4>q,x^). 


If linearly independent, let “ ^2 

[N( 4 >^,x-)] = 0 

(13) 

[N(x^,X 2 >] = 0 , 

if they are dependent, continue on with ••• ••• until 

determining a vector which is linearly independent. In this manner, 

it is possible .;o generate a set of linearly independent vectors 
(x^, ... »x^)^ £ n- *■ i» from a finite set of vectors 

(C^, ... , €^), m ^ k, belonging to > i = 1,2 By construction 

there exist constants ^ “ 1 , 2 ,... , m = 1,2 k such that 

£ , = a . , X, + ... + a X, 
j jl 1 jk k 

for all j = 1 , 2 , . . . , 

Now the conditions (12), (13), ... can be summarized as 


[N(4>q,x^)] = 0 

[N(x^,x^)] = 0 

for i,j = 1 , 2 ,. ..,k , i < 3 . 


Equation (I4) along with [N(4>Q,(j)Q) ] = 0 consists of 


(k^ + 3k + 4) < I (n^ + n + 2) 



conditions for k < n - 1. 


Now N(x,x) can be written as 


N(x,x) = (N.x, ... ,N x)x 
i n 

where N^, ... ,N^ are symmetric n x n matrices, which yields 
12 

Y n (n + 1) elements which can be chosen arbitrarily to satisfy (13) . 
12 12 

Since ^ ^ ^ ^ + n + 2) , equation (13) represents a linear 

12 

system in the ^ ^ (n + 1) arbitrary elements of the ... ,N^ , 

it is clear that conditions can be imposed to assure the existence of 
the quadratic nonlinear term from (13). 

Remark . If conditions are imposed on the solvability condition (9) 
to insure that o(e) ^ 0, the vertical bifurcation does not occur. 
Further, using Sattinger's arguments on the exchange of stability, 
it is possible to prove that the bifurcating branch of solutions now 
obtained by the implicit function theorem is either asymptotically stable 
or unstable from the first variation equation of (1). Appealing to 
a theorem of iiartman [5], we can now conclude that the bifurcating branch 
of solutions is structurally stable for e > 0, e sufficiently small. 



^3 

References 

[1] M, A. Krasnosel' skii. Topological Methods in the Theory o£ Ro'alir.ear 
Integral Equations, Macmillan, New York, 1S64. 

[2] M. W, Eirsch, S. Smale, Differential Equations, Dynamical Systems 
and Linear Algebra, Academic Press, New York, 1974, 

[3] D. H. Satcinger, Topics in Stability and Bifurcation Theory, Springer 
Lecture Notes in Mathematics, 309, Berlin, 1973. 

[4] D. ,i. Sattinger, Six Lectures on the Transition to Instability, 

Nonlinear Problems in the Physical Sciences and Biology, I. Stakgold, ed., 
Lecture Notes in Mathematics, 322, Berlin, 1973. 

[5] P. Earcma-n, A Lemma in the Theory of Structural Stability of Differential 
Equatroas, Proc. Amer. Math. Soc. IT (1960), 610-620. 


y 



APPENDIX XIII • 


SCALAR n-th ORDER COMPARISON SYSTEMS 


R. W. Gunderson* 


Department of Mathematics, Utah State University, Logan, Utah 84321. This 
work was supported by the National Aeronautics and Space Administration 
under Grant NGR-45-002-016 . 



INTRODUCTION 


The theory of differential inequalities has found an important application 
in the theory of comparison systems (vector Liapunov functions) . The theory 
of composite systems, for example, has been developed extensively since first 
introduced by Bailey in [1] (see Siljak [2] for an excellent survey) . At 
the same time, there have been several papers dealing with the application of 
higher order trajectory derivatives in the analysis of dynamical systems [3], 
[^]> [5], and [6]. In this note we shall discuss some recent results on the 
application of differential Inequalities to comparison theorems involving 
higher order trajectory derivatives. 

NOTATION AND DEFINITIONS 

Consider the system of first order differential equations 

X = f(t,x) (1) 

and the m-th order comparison system 

u^'"^(t) = u)(t,u,u’ , . . . .u^™ ^^) (2) 

where x, f belong to r” and t, u are scalars. Assume f continuous on 

D^ = {(t,x) 0 t T < X < r} and the right side of (2) continuous 

on [0,T] X r"*. a solution of (1) satisfying the initial condition x at 

o 

t^ will be denoted by x(t,t^,x^) and a solution of (2) satisfying 

u^ (t^) = Uj (j = l,2,...,m-l) will be denoted by u(t,t^,U^). 

Definition . The scalar function g(x), x belonging to r”, will be said 

to be of type W* on a set S r’^ if g(a) ^ g(b) for any a, b in S such that 

a = b , a. < b, (i=l,2, . . . ,n-l) . 
n n i— i 

COMPARISON THEOREMS 
The following lemma was given in [5] : 

Lemma. Let v;D^ — ^ R and let v belong to c”*, f belong to c"* ^ on D^. Let (o 



of 


equation (2) be of type W* in S for each t, where S = {t, 
V™ ^ (t,x)} and 


v^^^(t,x) = 


9v 

3t 


(j-1) 


9v 

9x 


(j-1) 


f (t,x) 


v(t,x) ,v' (t,x) , 
(3) 


Suppose 

v^"’^(t,x) £ o)(t,v,v’ , . . . (4) 

for (t,x) belonging to and set v^(0,x^) = u^ . Let J denote the maximal 
Interval of existence of the right maximal solution u (t,0,U ). Then 

v^^\t,x(t,0,x^)) £ u^^\t,0,U^) (5) 

for each t belonging to J [0,T] and j=l,2, . . . ,m-l. 

Remark . Since the inequality (4) implies the system of first order 
inequalities 


V — V 
1 2 


^2 “ "^3 


( 6 ) 


V < Cj(t,v ,v , . . . ,v ) 
n “ 12 m 


it is natural to view higher order derivatives as a means of obtaining 
vector Liapunov functions. (The fact that m is of type W* assures that 
the right side of (6) satisfies the corresponding monotone property for 
vector functions). There are several practical difficulties, however, to 
this approach of constructing vector Liapunov functions. For example, one 
such difficulty is the requirement that the right side of (4) be of type 
W*. Since the characteristic polynomial of a linear n-th order comparison 


equation must then be of the form 


J. *^-l 

s+a,s -a_s 
n-1 n-2 


n-2 


- . . . - a.s - a = 0 (a. > 0) 

1 o 1 — 


to be of type W*, it follows that even in the linear case, the comparison 
equation cannot be stable, except when n=l. 

In order to improve the applicability of scalar n-th order comparison 


theorems, we introduce the following special case of a more general theorem 
whose proof is given in [7], 



Theorem 1. Assume 




(n) f. , (n-1), 

X =U)(t,X,X ,...,X 


(7) 


where ^ is the right side of (2), holds for t belonging to some interval I. 

Suppose the solution z(t,t,E) of the adjoint equation to (7) is non-negative 

on the interval t , t for each t > t ; t , t belonging to I, and where 
o — o o 

E = (1,0,..., 0). Then (7) and (2) together imply 

u(t) < x(t,t ,U ) 

— o o 

for all t ^ t belonging to I. 

The following is an example of a second order system which is not of 
type W* but which satisfies the conditions of the theorem: 

Example . .The second order system 

x" + 3x' + 2x = 0 


has the adjoint equation 


z" - 3x' + 2z = 0 


which has the solution 


satisfying the conditions 


/.s (t-t ) (t-t ). 

z(t) = e o (2 - e o ) 


z(t^) = 1 

z ' (t ) = 0 
o 

Clearly z(t) > 0 for all t < t^, so that the conclusion of the theorem 
holds on any interval I. 

In fact, it is not difficult to show that the conditions of the theorem 
are satisfied by any second order system with real, unequal, roots. 

Using the above theorem it is then possible to prove comparison theorems 
of the following type [7]: 

T 

Theorem 2 . Let v(x) = x Hx be a positive definite quadratic form and 

suppose the trajectory derivatives of v formed relative to (1) satisfy 

(m-1) 


(m) , 

V + a , u 
m-1 


+ . . . + a V < 0 
o — 
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